
*****************************************************************************************************
******Analysis 2: Liberal vs. Conservative

******Effects of Information Provision (Without considering misperception): Panel A of Figure 4

use analysis.dta, clear

drop if ideo5 == 6

reg Q19_concernforyou treatment NEP if ideo3  == 1
estimates store concern_self_L

reg Q19_concernforyou treatment NEP if ideo3  == 3
estimates store concern_self_C


reg Q20_concernfornation treatment NEP if ideo3  == 1
estimates store concern_nation_L

reg Q20_concernfornation treatment NEP if ideo3  == 3
estimates store concern_nation_C


reg Q23_1_normgov treatment NEP if ideo3  == 1
estimates store norm_gov_L

reg Q23_1_normgov treatment NEP if ideo3  == 3
estimates store norm_gov_C


reg Q23_2_normperson treatment NEP if ideo3  == 1
estimates store norm_person_L

reg Q23_2_normperson treatment NEP if ideo3  == 3
estimates store norm_person_C


reg consumer_behav treatment NEP if ideo3  == 1
estimates store consumer_L

reg consumer_behav treatment NEP if ideo3  == 3
estimates store consumer_C


reg willing_pay treatment NEP if ideo3  == 1
estimates store willingness_L

reg willing_pay treatment NEP if ideo3  == 3
estimates store willingness_C

reg asso_resp treatment NEP if ideo3  == 1
estimates store group_L

reg asso_resp treatment NEP if ideo3  == 3
estimates store group_C


reg poli_behav treatment NEP if ideo3  == 1
estimates store political_L

reg poli_behav treatment NEP if ideo3  == 3
estimates store political_C   

#delimit ;
coefplot 
    (concern_self_L, keep(treatment) aseq(Concern for self and family) \
     concern_nation_L, keep(treatment) aseq(Concern for the country) \
     norm_gov_L, keep(treatment) aseq(Government should do more) \
     norm_person_L, keep(treatment) aseq(Personal obligation) \
     consumer_L, keep(treatment) aseq(Consumer behaviors) \
     willingness_L, keep(treatment) aseq(Policy support/Willingness to pay) \
     group_L, keep(treatment) aseq(Group participation and contribution) \
     political_L, keep(treatment) aseq(Political activities) \, 
     mcolor(gs1) msymbol(O) msize(vsmall) 
     ciopts(color(gs1 gs1) lcolor(*.5 *.5)) offset(0.10)) 
 
     
    (concern_self_C, keep(treatment) aseq(Concern for self and family) \
     concern_nation_C, keep(treatment) aseq(Concern for the country) \
     norm_gov_C, keep(treatment) aseq(Government should do more) \
     norm_person_C, keep(treatment) aseq(Personal obligation) \
     consumer_C, keep(treatment) aseq(Consumer behaviors) \
     willingness_C, keep(treatment) aseq(Policy support/Willingness to pay) \
     group_C, keep(treatment) aseq(Group participation and contribution) \
     political_C, keep(treatment) aseq(Political activities) \, 
     mcolor(gs13) msymbol(S) msize(vsmall) 
     ciopts(color(gs13 gs13) lcolor(*.5 *.5)) offset(-0.10)) 
     
     ,swapnames 
     levels(95 90) xline(0, lwidth(thin) lpattern(dash) lcolor(black)) 
     coeflabels(,notick labgap(2)) ylabel(,labsize(small)) plotlabels("Liberal" "Conservative") 
     xtitle("Marginal Effect of Treatment", size(small)) 
     headings("Concern for self and family" = "{bf:Issue-specific Environmental Concern}" "Government should do more" = "{bf:Norms/Attitudes}" "Consumer behaviors" = "{bf:Preferences and Behavioral Intentions}", nogap labgap(2)) 
     graphregion(fcolor(white)) grid(none) legend(size(small) rows(1))
;
#delimit cr


****** Effects of Information Provision by misperception: Panels B-D of Figure 4

use analysis.dta, clear

drop if ideo5 == 6
keep if misbelief > -80 & misbelief < 40

reg Q19_concernforyou i.treatment c.misbelief i.treatment#c.misbelief NEP if ideo3 == 1
estimates store concern_self_L_reg
margins, dydx(treatment) at(misbelief=(-80(5)40)) post
estimates store concern_self_L

reg Q19_concernforyou i.treatment c.misbelief i.treatment#c.misbelief NEP if ideo3 == 3
estimates store concern_self_C_reg
margins, dydx(treatment) at(misbelief=(-80(5)40)) post
estimates store concern_self_C



reg Q20_concernfornation i.treatment c.misbelief i.treatment#c.misbelief NEP if ideo3 == 1
estimates store concern_nation_L_reg
margins, dydx(treatment) at(misbelief=(-80(5)40)) post
estimates store concern_nation_L

reg Q20_concernfornation i.treatment c.misbelief i.treatment#c.misbelief NEP if ideo3 == 3
estimates store concern_nation_C_reg
margins, dydx(treatment) at(misbelief=(-80(5)40)) post
estimates store concern_nation_C

 

reg Q23_1_normgov i.treatment c.misbelief i.treatment#c.misbelief NEP if ideo3 == 1
estimates store norm_gov_L_reg
margins, dydx(treatment) at(misbelief=(-80(5)40)) post
estimates store norm_gov_L

reg Q23_1_normgov i.treatment c.misbelief i.treatment#c.misbelief NEP if ideo3 == 3
estimates store norm_gov_C_reg
margins, dydx(treatment) at(misbelief=(-80(5)40)) post
estimates store norm_gov_C


reg Q23_2_normperson i.treatment c.misbelief i.treatment#c.misbelief NEP if ideo3 == 1
estimates store norm_person_L_reg
margins, dydx(treatment) at(misbelief=(-80(5)40)) post
estimates store norm_person_L

reg Q23_2_normperson i.treatment c.misbelief i.treatment#c.misbelief NEP if ideo3 == 3
estimates store norm_person_C_reg
margins, dydx(treatment) at(misbelief=(-80(5)40)) post
estimates store norm_person_C


reg consumer_behav i.treatment c.misbelief i.treatment#c.misbelief NEP if ideo3 == 1
estimates store consumer_L_reg
margins, dydx(treatment) at(misbelief=(-80(5)40)) post
estimates store consumer_L

reg consumer_behav i.treatment c.misbelief i.treatment#c.misbelief NEP if ideo3 == 3
estimates store consumer_C_reg
margins, dydx(treatment) at(misbelief=(-80(5)40)) post
estimates store consumer_C


reg willing_pay i.treatment c.misbelief i.treatment#c.misbelief NEP if ideo3 == 1
estimates store willingness_L_reg
margins, dydx(treatment) at(misbelief=(-80(5)40)) post
estimates store willingness_L

reg willing_pay i.treatment c.misbelief i.treatment#c.misbelief NEP if ideo3 == 3
estimates store willingness_C_reg
margins, dydx(treatment) at(misbelief=(-80(5)40)) post
estimates store willingness_C

reg asso_resp i.treatment c.misbelief i.treatment#c.misbelief NEP if ideo3 == 1
estimates store group_L_reg
margins, dydx(treatment) at(misbelief=(-80(5)40)) post
estimates store group_L

reg asso_resp i.treatment c.misbelief i.treatment#c.misbelief NEP if ideo3 == 3
estimates store group_C_reg
margins, dydx(treatment) at(misbelief=(-80(5)40)) post
estimates store group_C


reg poli_behav i.treatment c.misbelief i.treatment#c.misbelief NEP if ideo3 == 1
estimates store political_L_reg
margins, dydx(treatment) at(misbelief=(-80(5)40)) post
estimates store political_L

reg poli_behav i.treatment c.misbelief i.treatment#c.misbelief NEP if ideo3 == 3
estimates store political_C_reg 
margins, dydx(treatment) at(misbelief=(-80(5)40)) post
estimates store political_C   

#delimit ;   
coefplot
    (concern_self_L, at label(Liberal) recast(line) lcolor(gs1) lwidth(*1.5) ciopts(recast(rline) color(gs1) lpattern(shortdash_dot) lcolor(%85)))
    (concern_self_C, at label(Conservative) recast(line) lcolor(gs13) lwidth(*1.5) ciopts(recast(rline) color(gs13) lpattern(dash) lcolor(%85))), 
     bylabel(Concern for Self and Family)
    || 
    (concern_nation_L, at recast(line) lcolor(gs1) lwidth(*1.5) ciopts(recast(rline) color(gs1) lpattern(shortdash_dot) lcolor(%85)))
    (concern_nation_C, at recast(line) lcolor(gs13) lwidth(*1.5) ciopts(recast(rline) color(gs13) lpattern(dash) lcolor(%85))),
    bylabel(Concern for the Country)
    ||, 
     byopts(title("Issue-specific Environmental Concern", size(medium)))

     ytitle("Marginal Effect of Treatment", size(small)) 
     xtitle("Misperception (Perceived Ranking of Risk - Actual Ranking of Risk)", size(small))
     yline(0, lwidth(thin) lpattern(dash) lcolor(black))
     graphregion(fcolor(white))
     legend(size(small) rows(1))
;
#delimit cr


#delimit ;   
coefplot
    (norm_gov_L, at label(Liberal) recast(line) lcolor(gs1) lwidth(*1.5) ciopts(recast(rline) color(gs1) lpattern(shortdash_dot) lcolor(%85)))
    (norm_gov_C, at label(Conservative) recast(line) lcolor(gs13) lwidth(*1.5) ciopts(recast(rline) color(gs13) lpattern(dash) lcolor(%85))), 
     bylabel(Government Should Do More) legend(size(vsmall) rows(1))
    || 
    (norm_person_L, at recast(line) lcolor(gs1) lwidth(*1.5) ciopts(recast(rline) color(gs1) lpattern(shortdash_dot) lcolor(%85)))
    (norm_person_C, at recast(line) lcolor(gs13) lwidth(*1.5) ciopts(recast(rline) color(gs13) lpattern(dash) lcolor(%85))),
    bylabel(Personal Obligation)
    ||, 
     byopts(title("Norms/Attitudes", size(medium)) legend(size(vsmall) rows(1)))

     ytitle("Marginal Effect of Treatment", size(small)) 
     xtitle("Misperception (Perceived Ranking of Risk - Actual Ranking of Risk)", size(small))
     yline(0, lwidth(thin) lpattern(dash) lcolor(black))
     graphregion(fcolor(white))
     legend(size(small) rows(1))
;
#delimit cr


#delimit ;   
coefplot
    (consumer_L, at label(Liberal) recast(line) lcolor(gs1) lwidth(*1.5) ciopts(recast(rline) color(gs1) lpattern(shortdash_dot) lcolor(%85)))
    (consumer_C, at label(Conservative) recast(line) lcolor(gs13) lwidth(*1.5) ciopts(recast(rline) color(gs13) lpattern(dash) lcolor(%85))), 
     bylabel(Consumer Behaviors) 
    || 
    (willingness_L, at recast(line) lcolor(gs1) lwidth(*1.5) ciopts(recast(rline) color(gs1) lpattern(shortdash_dot) lcolor(%85)))
    (willingness_C, at recast(line) lcolor(gs13) lwidth(*1.5) ciopts(recast(rline) color(gs13) lpattern(dash) lcolor(%85))),
    bylabel(Policy Support/Willingness to Pay)
    ||
    (group_L, at recast(line) lcolor(gs1) lwidth(*1.5) ciopts(recast(rline) color(gs1) lpattern(shortdash_dot) lcolor(%85)))
    (group_C, at recast(line) lcolor(gs13) lwidth(*1.5) ciopts(recast(rline) color(gs13) lpattern(dash) lcolor(%85))),
    bylabel(Group Participation and Contribution)
    ||
    (political_L, at recast(line) lcolor(gs1) lwidth(*1.5) ciopts(recast(rline) color(gs1) lpattern(shortdash_dot) lcolor(%85)))
    (political_C, at recast(line) lcolor(gs13) lwidth(*1.5) ciopts(recast(rline) color(gs13) lpattern(dash) lcolor(%85))),
    bylabel(Political Activities)
    ||, 
     byopts(title("Preferences and Behavioral Intentions", size(medium)) legend(size(vsmall) rows(1)))

     ytitle("Marginal Effect of Treatment", size(small)) 
     xtitle("Misperception (Perceived Ranking of Risk - Actual Ranking of Risk)", size(small))
     yline(0, lwidth(thin) lpattern(dash) lcolor(black))
     legend(size(small) rows(1))
     graphregion(fcolor(white))
;
#delimit cr